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CHAOTIC  DYNAMICS  OF  ELASTIC-PLASTIC  BEAMS 
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During  the  last  decades  great  attention  has  been  turned  to  chaotic  vibrations  of  elastic 
beams,  but  not  much  has  been  done  in  the  case  of  elastic -plastic  deformations.  Evidently 
the  first  paper  in  this  field  belongs  to  Symonds  and  Yu  (1985),  who  considered  the 
following  problem.  A fixed  ended  beam  is  subjected  to  short  intensive  pulse  of  transverse 
loading  that  produces  plastic  deformation.  Since  the  ends  of  the  beam  are  fixed  membrane 
forces  must  be  taken  into  account.  Solving  the  equations  of  motion  Symonds  and  Yu 
found  that  permanent  deflection  may  be  in  direction  opposite  to  the  load.  This 
phenomenon  was  investigated  in  several  papers  by  Symonds  and  his  collaborators.  It 
turned  out  that  the  permanent  defflection  is  very  sensitive  to  small  changes  of  load. 
Fractality  and  self-similarity  which  are  characteristic  to  chaotic  processes,  were 
demonstrated.  For  the  similarity  dimension  the  value  0,78  and  for  the  correlation  fractal 
dimension  ~1,44  were  obtained. 

In  most  papers  of  Symonds  and  his  collaborators  the  real  beam  was  replaced  by  a 
Shanley-type  model  with  one  or  two  degrees  of  freedom.  In  some  papers  (see  Symonds 
and  Qian,  1996)  the  Galerkin  method  for  2DoF  was  used;  some  simplifying  assumptions 
(sandwich  beam,  perfectly  plastic  material,  disregard  of  axial  inertia  forces)  were  made. 

The  aim  of  this  paper  is  to  present  a method  of  solution,  which  is  based  on  the 
Galerkin  technique  and  is  applicable  for  beams  with  an  arbitrary  number  of  DoF.  Basic 
elements  of  this  method  were  worked  out  by  Lepik  (1994,  1995). 

To  be  brief  we  shall  consider  beams  with  a rectangular  cross-section,  B,h  and  L are 
the  width,  thickness  and  length  of  the  beam,  respectively.  The  equations  of  motion  are 
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where  p is  density,  u - axial  displacement,  w - deflection.  Axial  force  T and  bending 

moment  M we  shall  calculate  from  the  formulae 
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We  shall  assume  that  the  beam  material  has  linear  strain-hardening;  elastic  unloading 
and  secondary  plastic  loading  effects  are  also  taken  into  account. 

To  the  equations  (1)  we  shall  apply  the  Galerkin  procedure 
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In  the  case  of  a beam  with  simply  supported  ends  we  shall  seek  the  solution  in  the 
form 
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where  e(0)  = e0  and  the  coefficients  e0,ak,fk  are  subject  to  variation.  The  integrals  (2) 
will  be  calculated  numerically;  replacing  these  results  into  equations  (3)  we  find  the 
second  derivatives  e0 ,ak,fk.  The  quantities  e0 ,ak,fk  are  evaluated  according  to  the 
method  of  central  finite  differences  (for  getting  a stable  solution  the  time  increment  At 
must  be  sufficiently  small). 

In  order  to  demonstrate  the  efficiency  of  the  proposed  algorithm  several  numerical 
examples  are  examined.  With  the  purpose  to  establish  chaotic  effects  the  deflection 
history,  phase  and  power  spectrum  diagrams  are  put  together.  It  turned  out  that  chaotic 
motion  of  the  beam  may  exist,  especially  in  the  initial  phase  of  motion;  as  to  the  long 
term  motion  then  it  translates  to  periodic  vibrations  of  smaller  amplitude. 
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